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Abstract. The presence of macroscopic phase separation into superconducting and 
magnetic phases in LaFeAsOi_2;Fa; and CaFei_a;Coa;AsF is demonstrated by muon 
spin rotation (/iSR) measurement across their phase boundaries {x = 0.06 for 
LaFeAsOi-jjE^; and x = 0.075-0.15 for CaFei_3;Coa;AsF). In LaFeAsO0.94F0.06, both 
magnetism and superconductivity develop simultaneously below a common critical 
temperature, T^i ~ Tc — 18 K, where the magnetism is characterized by strong 
randomness. A similar, but more distinct segregation of these two phases is observed 
in CaFei-xCoxAsF, where the magnetic phase retains T,„ as close to that of the parent 
compound (Tc <C Tm — 80-120 K) and the superconducting volume fraction is mostly 
proportional to the Co content x. The close relationship between magnetism and 
superconductivity is discussed based on these experimental observations. Concerning 
superconducting phase, an assessment is made on the anisotropy of order parameter in 
the superconducting state of LaFeAsOi-ajF^^, CaFei-xCo^AsF, and Bai_2;K2:Fe2As2 
(x = 0.4) based on the temperature dependence of superfluid density [ns{T)] measured 
by /iSR. The gap parameter, 2A/kBTc^ determined from ns{T) exhibits a tendency 
that hole-doped pnictides (Bai_2:K^Fe2As2) is much greater than those in electron- 
doped ones (LaFeAsOi-xFx, CaFci-xCoajAsF), suggesting difference in the couphng 
to bosons mediating the Cooper pairs between relevant d electron bands. 
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1. Introduction 

The recent discovery of tlie iron pnictide superconductor LaFeAsOi_a;Fa; (LFAO-F) over 
a fluorine concentration of 0.05 < a; < 0.2 with the maximal critical temperature (T^) of 
26 K [l] and the following revelation of much increased Tc upon the substitution of La for 
other rare-earth elements (Ce, Pr, Nd, Sm,... leading to a maxium Tc of 55 K [2, [3l H]) 
or the application of pressure for LFAO-F (~43 K [5]) have triggered broad interest in 
the mechanism yielding a relatively high Tc in this new class of compounds. They have 
a layered structure like high-Tc cuprates, where the dopant and conducting layers are 
so separated that carriers {electrons, in this case) introduced by the substitution of 
with F~ in the La202 layers move within the layers consisting of strongly bonded Fe 
and As atoms. Moreover, very recent developments demonstrate increasing variety in 
the methods of electron doping such as oxygen depletion [HI [7] or Co substitution for 
Fe [H [9l [TOl [11]. A similar situation is presumed for the ternary compound y4Fe2As2 
{A = Ba, Sr, Ca), where holes are introduced by the substitution of A^'^ with 
ions {B = Na, K, Cs) [121 [131 [IH [IS]- They exhibit another qualitative similarity to 
cuprates in that superconductivity occurs upon carrier doping of pristine compounds 
that exhibit magnetism [l6l[I7l[l8l[T9l[2ni[2ll[22l[23]. Recent results of the muon spin 
rotation/relaxation (/xSR) experiment on a variety of iron pnictide superconductors 
showed that the superfluid density rig may fall on the empirical line on the rig vs Tc 
diagram observed for the underdoped cuprates [TH [2l], from which possibility of the 
common mechanism of superconductivity is argued between oxypnictides and cuprates. 

The iron pnictides exhibit an interesting similarity with cuprates that the variation 
of Tc against doping is "bell-shaped" in hole-doped compounds [TH [13] while Tc does 
not vary much with x in electron-doped case [H [25]. However, recent investigations 
in electron-doped (n type) cuprates strongly suggest that such an electron-hole 
"asymmetry" is a manifestation of difference in the fundamental properties of underlying 
electronic states between these two cases, where the n type cuprates are much more like 
normal Fermi liquids rather than doped Mott insulators [26]. This might be readily 
illustrated by pointing out that, given all the doped carriers participate in the Cooper 
pairs (as suggested by experiment), the insensitivity of Tc against the variation of 
(oc x) observed in LFAO-F [1] cannot be reconciled with the above-mentioned empirical 
linear relation, while it is reasonably understood from the conventional BCS theory 
where condensation energy is predicted to be independent of carrier concentration. More 
interestingly, the very recent revelation of superconductivity upon Co substitution for 
Fe in LaFeAsO and other iron pnictides (where the Co atoms serves as electron donors) 
brings out a sheer contrast between these two classes of materials in terms of response 
to the substitution of transition metal ions as well as the tolerance of superconductivity 
to the distortion of conducting layers [8], [9l [TDl [H] . 

The close relationship of magnetism and superconductivity suggests that a detailed 
investigation of how these two phases coexist (and compete) near the phase boundary 
will provide important clues to elucidating the paring mechanism. Among various 
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techniques, /xSR has a great advantage in that it can be apphed in systems consisting 
of spatially inhomogeneous multiple phases, providing information on respective phases 
according to their fractional yield. Our /iSR measurement in a LFAO-F sample with 
X = 0.06 (Tc ~ 18 K) reveals that these two phases indeed coexist in the form of 
macroscopic phase separation, and more interestingly, that a spin glass-like magnetic 
phase develops in conjunction with superconductivity in the paramagnetic phase |27j . 
This accordance strongly suggests a common origin of the electronic correlation between 
these two competing phases. On the other hand, /iSR study on CaFei_j,Coa;AsF 
(CFCAF, a variation of LFAO-F with trivalent cation and oxygen respectively replaced 
with divalent alkali metal and fluorine, and the carrier doping is attained by substituting 
Co for Fe) reveals a unique character of the Fe2As2 layers that the superconducting 
state is realized over a vicinity of Co atoms, as inferred from the observation that the 
superconducting volume fraction is nearly proportional to the Co concentration |28j . 
The rest of the CFCAF specimen remains magnetic (strongly modulated spin density 
wave), thus indicating that superconductivity coexists with magnetism again in a form 
of phase separation. 

As already mentioned, /iSR can provide information selectively from the 
superconducting parts of samples even in the situation that they coexist with 
magnetism. We examine the temperature dependence of in LFAO-F, CFCAF, 
and in Bai_a;Ka;Fe2As2 (BKFA) |29], and discuss the degree of anisotropy in their 
superconducting order parameters and the strength of coupling to bosons that mediate 
the Cooper pairs. 

2. Experiment 

Although the oxypnictides i^FeAsOi-^^Fa, with rare-earth substitution {R = Nd, Sm, 
etc.) exhibits higher than that of LFAO-F, strong random magnetic fields from 
rare-earth ions preclude a detailed study of the ground state using sensitive magnetic 
probes hke //SR. In this regard, the original LFAO-F (as well as other two compounds 
in this review) has a major advantage that we can attribute the origin of magnetism, 
if at all detected by /iSR, to that of the Fe2As2 layers without ambiguity. The target 
concentration has been chosen near the phase boundary, x = 0.06 for LaFeAsOi-ajF^. 
and 0.075-0.15 for CaFei„i.Co2:AsF, where a polycrystalline sample has been synthesized 
by solid state reaction. The detailed procedure for sample preparation is described in 
an earlier report [U [11], EH] • For Bai_xKa;Fe2As2, polycrystalline samples with nominal 
compositions of x = 0.1 and 0.4 were prepared by a solid state reaction using starting 
materials having the highest purity available, where details of the preparation process 
are described in an earlier report The obtained LFAO-F sample was confirmed to 
be mostly of single phase using X-ray diffraction method (see figure [1]). Of two possible 
impurity phases, namely LaOF and FeAs, only the latter exhibits a magnetic (helical) 
order with T/v ^ 77 K [301 [H] ■ For the CFCAF samples, it was inferred from detailed X- 
ray diffraction analysis that nearly 11% of the sample was crystalized in Ca(Fe,Co)2As2 
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for the sample with x = 0.15, while no such phase was observed for other samples. In 
addition, a small fraction (2-4%) of CaF2 (fluorite) was found to exist as impurity in all 
samples. Since muons in fluorite are depolarized immediately upon implantation, their 
contribution should be negligible. The influence of Ca(Fe,Co)2As2 will be discussed in 
the following section. Finally, all the observed diffraction peaks in the Bai_2;K2;Fe2As2 
sample with x = 0.4 was perfectly reproduced by those of single phase compound [T2]. 
while a minor unknown impurity phase (~2%) was observed for x = 0.1. 

Conventional /iSR measurement was performed using the LAMPF spectrometer 
installed on the M20 beamline of TRIUMF, Canada. During the measurement under a 
zero field (ZF), residual magnetic field at the sample position was reduced below 10~^ T 
with the initial muon spin direction parallel to the muon beam direction [Pfi{0) \\ z]. 
For longitudinal field (LF) measurement, a magnetic field was applied parallel to -P/x(0). 
Time-dependent muon polarization [Gz{t) = z ■ P/iit)] was monitored by measuring 
decay-positron asymmetry along the z-axis. Transverse field (TF) condition was realized 
by rotating the initial muon polarization so that -P/^(0) || x, where the asymmetry was 
monitored along the x-axis to obtain Gx{t) = x ■ P^{t). All the measurements under 
a magnetic field were made by cooling the sample to the target temperature after the 
field equilibrated. 

3. Result 

3.1. LaFeAsOi^^F^ 

The quality of samples measured by /xSR has been examined by looking into X- 
ray diffraction and magnetic susceptibility data. As shown in figure [H magnetic 
susceptibility exhibits no trace of FeAs phase or local magnetic impurities except below 
~ 50 K where a small upturn is observed. The susceptibility at a lower field [shown in 
figure El^a)] provides evidence of bulk superconductivity with Tc ~ 18 K from the onset 
of diamagnetism. 

ZF-/iSR is the most sensitive technique for examining magnetism in any form, where 
the development of local magnetic moments leads to either the spontaneous oscillation 
(for long-range order) or exponential damping (inhomogeneous local magnetism) of 
Gz{t). Figure [2] shows examples of ZF-/iSR time spectra collected at 2 K and 30 K. 
The spectrum at 30 K (> T^) exhibits a Gaussian-like depolarization due to weak 
random local fields from nuclear magnetic moments, indicating that the entire sample 
is in paramagnetic state. Meanwhile, the spectrum at 2 K is split into two components, 
one that exhibits a fast depolarization and the other that remains to show Gaussian- 
like relaxation. This indicates that there is a finite fraction of implanted muons that 
sense hyperfine fields from local electronic moments. The absence of oscillatory signal 
implies that the hyperfine field is highly inhomogeneous, so that the local magnetism 
is characterized by strong randomness or spin fluctuation. The fractional yield of this 
component is as large as 25% (see below), which is hardly attributed to impurity and 
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Figure 1. (Color online) X-ray diffraction (top) and Magnetic susceptibility data 
obtained for the sample of LaFeAsO0.94F0.06 rneasured by /iSR measurement. Inset 
shows a reduced view of the region below 35 K. 

therefore implies that the sample exhibits a macroscopic phase separation into two 
phases. 

The magnitude of the hyperfine field and that of spin fiuctuation are evaluated 
by observing the response of the /iSR spectrum to a longitudinal field (LF). It is 
shown in figure [2] that the slow depolarization of the paramagnetic component is 
quenched by applying a weak magnetic field (LF=5 mT), which is perfectly in line 
with the suppression of static nuclear dipolar fields (< 10° mT). Meanwhile, the faster 
depolarization (seen for < t < 1 /is) due to the magnetic phase is recovered only 
gradually over a field range of 10^'"^ mT, and there still remains a slow depolarization 
even at the highest field of 60 mT. This residual depolarization under LF is a clear sign 
that local spins are slowly fiuctuating, leading to the spin- lattice relaxation of P^{t). 
Such quasi-two-step depolarization is also observed in dilute spin-glass systems. [32] 
which is understood as a distribution of spin correlation time. A detailed analysis 
is made considering that these two components coming from the magnetic phase (see 
below). 

Under a transverse field, implanted muons experience an inhomogeneity of the 
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Figure 2. (Color online) /iSR time spectra observed in LaFeAsOi_j;Fa; with x = 0.06 
at 2 K under a longitudinal field (LF), a zero field (ZF), and a transverse field (TF), 
and that under ZF at 30 K. The spectrum under TF is plotted on a rotating reference 
frame of 6.78 MHz to extract the envelop function. 
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0.165(9) 


0.71(5) 


1.7(2) 
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0.081(4) 


3.9(3) 


4(1) 



Table 1. Physical parameters obtained from LF-^SR spectra at 2 K by analysis using 
eq. O. 

field [i?2(r)] due to flux line lattice formation below Tc that leads to depolarization, in 
addition to those observed under a zero field. The TF-/iSR time spectrum in figure [2] 
(envelop part of the oscillation) obtained under a field of 50 mT exhibits complete 
depolarization at 2 K, indicating that the entire volume of the paramagnetic phase falls 
into the superconducting state. The rapidly depolarizing component observed under ZF 
is also visible (although the coarse binning of the spectra for extracting the envelop part 
makes it slightly obscure), indicating that the corresponding part of the sample remains 
magnetic. 

Considering the presence of the magnetic phase besides the paramagnetic 
(=superconducting below Tc) phase, we take special precaution to analyze both TF 
and ZF/LF /iSR spectra in a consistent manner. For the determination of physical 
parameters describing the behavior of signals from the magnetic phase, we first analyze 
ZF/LF spectra at 2 K using the x-square minimization fit with the relaxation function 

3 

G,{t) = [wi + J2 exp i-Kt)] ■ Gkt{S^ : t), (1) 

i=2 

where Gkt(5n : t) is the Kubo-Toyabe relaxation function for describing the Gaussian 
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damping due to random local fields from nuclear moments (with 6-^ being the linewidth) 
[33] , Wi is the fractional yield for the paramagnetic phase, W2 and 1^3 are those for the 
magnetic phase (Y^Wi = 1) with A2,3 being the corresponding relaxation rate described 
by the Redfield relation 

^^ = ^^ (^ = 2,3), (2) 

where iv^ = 7^ifLF, 7/i is the muon gyromagnetic ratio (= 2n x 135.53 MHz/T), i^Lp 
is the longitudinal field, 62 and ^3 are the mean values of the hyperfine fields being 
exerted on muons from local electronic moments, and i>2 and are the fluctuation 
rates of the hyperfine fields. The solid curves in figure [2] show the result of analysis 
where all the spectra at different fields (only ZF and LF=5mT are shown here) are 
fitted simultaneously using eqs. ([T]) and ([2]) with common parameter values (except 
that is fixed to the respective value for each spectrum), which show excellent agreement 
with all the spectra. The deduced parameters are shown in table [H Although the 
depolarization in the magnetic phase is approximately represented by two components 
with different hyperfine couplings (Si), the fluctuation rates (z/j) are close to each other 
(10'' s~^ at 2 K), suggesting that the randomness is primarily due to the distribution 
in the size of local moments (or in their distances to muons). Since no impurity phase 
with a fraction as large as 25% is detected by X-ray diffraction analysis, it is concluded 
that this magnetic phase is intrinsic. 

In the analysis of temperature-dependent TP spectra, we used the relaxation 
function 

G,{t) = [wie-^''' cos (27r/,t + (P) + {w2+W3)e-^^^' cos (27r/^t + (f))]e-^^'^'\ (3) 

where Wi and 5n are fixed to the values obtained by analyzing ZF/LF-/iSR spectra. 
The first component in the above equation represents the contribution of flux line 
lattice formation in the superconducting phase, where 6s corresponds to the linewidth 
(Ts = = 7/i((-B(r) — Bq)'^)^^'^ [with Bq being the mean B{r)], while the second 

term represents the relaxation in the magnetic phase. Here, the relaxation rate for the 
latter is represented by a single value Am (instead of A2,3), as it turns out that the two 
components observed under LP are hardly discernible in TF-/iSR spectra. [This does 
not affect the result of the analysis, because the amplitude is fixed to W2 + W3 so that A^ 
may represent a mean ~ {W2A2+W3A3) / {W2+W3).] The fit analysis using the above form 
indicates that all the spectra are perfectly reproduced while the partial asymmetry is 
fixed to the value determined from ZF-/iSR spectra. This strengthens the presumption 
that the paramagnetic phase becomes superconducting below T^. The result of analysis 
is summarized in figure [3l together with the result of dc magnetization measured in the 
sample from the same batch as that used for /iSR. 

It is interesting to note in figure [3](b) that, although the central frequency in the 
superconducting phase (/s) does not show much change below Tc ^ 18 K probably 
owing to a large magnetic penetration depth (it is indeed large, see below), that in the 
magnetic phase (/m) exhibits a clear shift in the negative direction below Tm ~ Tc. 
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Figure 3. (Color online) Temperature dependence of dc magnetic susceptibility 
measured at 1 mT (a), and that of physical parameters deduced by analyzing TF- 
liSR spectra in superconducting (/s, Ss) and magnetic (/m, Am) phases (b-c), and of 
CTs (= y/2Ss) proportional to superfluid density (d). Curves in (d) are fits by models 
described in the text. 



Table 2. Physical parameters deduced from curve fitting for the data in figure [3I^d). 
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0.73(6) 








5.19 
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The magnitude of the shift is as large as ~ 1% and thus is readily identified despite a 
relatively low external field of 50 mT. As shown in figure [3]^c), the relaxation rate in 
the magnetic phase (Am) also develops below Tc in accordance with the frequency shift, 
demonstrating that a spin-glass-like magnetism sets in below T^. Here, we note that the 
development of magnetic phase is already evident in the ZF/LF-/iSR spectra, and results 
are fully consistent with each other. The onset of superconductivity below is also 
confirmed by an increase of 6s, as observed in figure [3]^c). This remarkable accordance 
of onset temperature between magnetism and superconductivity strongly suggests that 
there is an intrinsic relationship between the superconducting and magnetic phases that 
leads to a common characteristic temperature. 

The temperature dependence of (Xs in figure [3](d) is compared with theoretical 
predictions for a variety of models with different order parameters. The weak-coupling 
BCS model (s-wave, single gap) apparently fails to reproduce the present data, as they 
exhibit a tendency to vary with temperature over the region T/T^ < 0.4. As a next 
possibility, we compare the data with a phenomenological double-gap model [Ml 135] . 

as(T) = a(0) - w6aiA,, T) - (1 - w)M^2, T), (4) 
6aiA, T) = ^l^°^ f{s, T).\l- f\e, T)]de, 

where Aj = 1 and 2) are the energy gap at T = 0, w is the relative weight for ^ = 1, 
k-Q is the Boltzmann constant, f{e,T) is the Fermi distribution function, and A(T) is 
the standard BCS gap energy. A fit using the above model shown by a dotted curve 
seems to exhibit reasonable agreement with data (the value of 2Ai/fcBTc in ref. [27] is 
incorrect and it should be replaced with that in table 2). These observations suggest 
that the superconducting order parameter is not described by a s-wave symmetry with 
single gap. When a power law, ag = cro[l — (T/Tc)^], is used in fitting the data, we 
obtain a curve shown by the broken line in figure [3](d) with an exponent /3 ~ 2. This 
seems to favor the recent theory of fully gapped anisotropic ±s-wave superconductivity 
that predicts more enhanced quasiparticle excitations associated with the smaller gap 
[36]. 

In the limit of extreme type II superconductors [i.e., A/^ ^ 1, where A is the 
effective London penetration depth and ^ = a/ $o/(2vrifc2) is the Ginzburg-Landau 
coherence length, $o is the flux quantum, and Hc2 is the upper critical field], Us is 
determined by A using the relation [37] 

= 2.74 X 10-^(1 - h) [l + 3.9(1 - hf] %\-', (5) 

where h = Hi[f/Hc2 and H^p is the magnitude of external field. While eq. ([5]) is a 
good approximation under limited conditions (as it introduces a relatively large cutoff 
~ \/2/^ in calculating the second moment, see ref. [38] for more detail), we would 
continue using it to keep the consistency with earlier analysis. From cxs extrapolated 
to T = and taking Hc2 ~ 50 T (ref. [39]), we obtain A=595(3) nm. Because of 
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Figure 4. (Color online) Temperature dependence of electrical resistivity in 
CaFei„^CO:rAsF with a; = 0, 0.075, and 0.15. 

the large anisotropy expected from the layered structure of this compound, A in the 
polycrystalline sample would be predominantly determined by in-plane penetration 
depth (Aab)- Using the formula A = l.SlAab for such a situation [IQ], we obtain 
Aab=454(2) nm. This value coincides with that expected from the aforementioned 
empirical relation between A'^^^ superconductors [121 [IS!- However, this may not be 
uniquely attributed to the superfluid density because A depends not only rig but also on 
the effective mass, ds oc A~^ = n^e^ /m*c^. 

Finally, we discuss the feature of the spin glass-like phase. Assuming that the local 
moments are those of off-stoichiometric iron atoms with a moment size close to that in 
the SDW phase (~ 0.36/iB [E]), the mean distance between muon and iron moments 
in the relevant phase is estimated to be ~ 0.5 nm from an average of 5j. Given the unit 
cell size (a = 0.403 nm, c = 0.874 nm [1]), this would mean that more than a quarter 
of iron atoms in the magnetic phase (i.e., ~ 7% of the entire sample) should serve as 
a source of local moments. It is unlikely that such a significant fraction of iron atoms 
remains as impurities in the present sample. 

It might also be noteworthy that there is an anomaly near Tin2 — 12 K in 
the susceptibility [the onset of ZFC/FC hysteresis in figure [1] and a steplike kink in 
figure [3]^a)]. This seems to be in accordance with the onset of a steeper increase in Am 
below Tni2, suggesting a change in magnetic correlation. 

3.2. CaFei^^Co^AsF 

/iSR measurements has been made on samples with x = 0, 0.075, and 0.15, where the 
pristine compound exhibits anomaly around 120 K in the resistivity while the latter two 
fall into superconducting state below Tc ~ 18 K and 21 K (defined as the midpoint of 
fall in the resistivity, see figure HI) |llj . The sample quality is examined by detailed X- 
ray diffraction analysis. As shown in figure O the peaks representing the Fe2As2 layers 
exhibit a progressive shift with Co content, which is consistent with the formation of 
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Figure 5. (Color online) X-ray diffraction pattern observed in CaFei-a^Coj^AsF with 
X = 0, 0.075, and 0.15. The [006] and [214] peaks exhibit a progressive shift with 
doping concentration, indicating that the majority of the samples consist of a uniform 
solid-solution. 



Table 3. Result of X-ray diffraction analysis on CaFei^xCoxAsF samples. 



X = 





X = 0.075 


X = 0.15 


CaFeAsF 


92.50% 


Ca(Fe,Co)AsF 98.60% 


Ca(Fe,Co)AsF 


85.74% 


CaFs 


2.24% 


CaFa 1.40% 


CaFs 


3.71% 


FeAs 


5.26% 




Ca(Fe,Co)2As2 


10.54% 



uniform solid-solution. Meanwhile, a slight broadening of the [006] peak may suggest 
the fluctuation of lattice constants due to randomness introduced by Co substitution. 
The impurity phases deduced from X-ray diffraction analysis is summarized in table [3l 
where it is noticeable that nearly 11% of the sample with x = 0.15 is crystalized in 
Ca(Fe,Co)2As2 (Cal22 phase), while no such phase was observed for other samples. It 
is known that CaFe2As2 is antiferromagnetic below T/v = 171 K, where the transition 
is associated with structural change to orthorhombic [31]. So far there seems to be no 
report on the electronic property upon Co substitution for this compound. We therefore 
presume that muons stopped in the Cal22 phase does not contribute to the /xSR signal 
below 171 K. 

It is inferred from ZF-/iSR measurement that the anomaly observed in the sample 
with X = corresponds to the occurrence of magnetic phase below Tm — 120 K. As 
shown in figure O the /iSR spectra below exhibit a spontaneous oscillation with 
a well defind frequency that approaches /m — 25 MHz with decreasing temperature. 
This is an indication that implanted muons sense a unique internal magnetic field 
-Bm = 27r/m/7^ ^ 0.18 T. The magnitude of -Bm is in good agreement with earlier 
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Figure 6. (Color online) ZF-/iSR time spectra observed in CaPeAsF (undoped) at 
140 K and 2 K. (inset) Frequency of spontaneous oscillation {ly = /m) vs temperature. 

yuSR results in i?FeAsO |IHl HI [23], where a commensurate spin density wave (SDW) 
with a reduced moment of ~0.25/xb at the iron sites is suggested [I^. It is also inferred 
from LF-/iSR spectra that the internal field is static within the time scale of fiSR (< 10^^ 
s). 

It has been reported that the Co doping is quite effective to suppress the anomaly 
in the resistivity at T^; it virtually disappears at x ~ 0.1 where the superconductivity 
seems to be close to its optimum as suggested from the maximal Tc ~ 22 K |llj . 
However, ZF-/iSR measurements in the samples with x = 0.075 and 0.15 indicates 
that the superconductivity does not develop uniformly over the specimen. As shown 
in figure [TJ the time spectra exhibit a character similar to those observed in LFAO-F 
{x = 0.06), namely, they consist of two components, one showing rapid depolarization 
and other showing slow Gaussian damping with the relative yield of the latter increasing 
progressively with x. A closer look into the earlier time range of the spectra obtained 
for X = 0.075 indicates that the rapid depolarization corresponds to a strongly damped 
oscillation with a frequency roughly equal to f^. This, together with the common 
onset temperature for magnetism (Tm — 120 K), confirms that the signal comes from 
the SDW phase with a strong modulation due to Co doping. Very recent neutron 
diffraction experiment suggests a similar situation for x = 0.06, although they only 
observe a volume-averaged signal jl3]. As inferred from TF-/iSR measurements (see 
below), the rest of the specimen exhibits superconductivity below T^. 

Considering these features, ZF-/iSR spectra are analyzed by the x-square 
minimization fit with the relaxation function 

G,{t) = [w, + W2GM] GktI^n : t), (6) 

Gm{t) = 1 + ^e-^-* cos(27r/^t + 0), (7) 
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Figure 7. (Color online) fiSK time spectra observed in CaFei-^CoxAsF [x = 0.075 
(a) and 0.15 (b)] at 2 K under a longitudinal field (LF, open circles), a zero field 
(ZF, triangles), and a transverse field (TF, squares), and that under ZF above Tm 
(triangles). The spectrum under TF is plotted on a rotating reference frame to extract 
the envelop function. 



where wi is the fractional yield for the nonmagnetic phase, and W2 is that for the SDW 
phase (^Wi = 1) with being the relaxation rate for the spontaneous oscillation. 
The first term in the second component represents the spatial average of cos^ with 9 
referring to the angle between direction of the initial muon polarization and that of 
the static local field at the muon site, which equals 1/3 in polycrystalline specimen 
under zero external field. (This term would be also subject to depolarization in the case 
of fluctuating local field, as observed in LFAO-F). The fractional yields of respective 
components are shown in figure [HI 

The superconducting properties are extracted from TF-/xSR spectra by fitting 
analysis in a manner similar to the case of LFAO-F, where we use a simplified form 
of eq. (IHl) with the factor W2 + replaced with W2, 

G,{t) = [wic-'"'" cos(27r/st + 0) + w;2G^(^)]e-5'5a^^ (8) 

It is confirmed that the spectra are perfectly reproduced while the partial asymmetry 
(wj) is fixed to the value determined from ZF-/iSR spectra, supporting the presumption 
that the paramagnetic phase becomes superconducting below Tc. The result of analysis 
is summarized in figure [9l 

Unlike the case of LFAO-F, the magnetic component develops at temperatures 
much higher than the superconducting transition temperature (Tm > Tc), although Tm 
is considerably reduced to ~80 K for x = 0.15 (see figure [9]). Meanwhile, the oscillation 
observed in figure [6] disappears in both of the Co-doped samples (figure [3, ZF, 2K), 
indicating that the magnetic order is strongly modulated. This observation provides 
evidence that the phase separation is not simply due to aggregation of cobalt atoms 
upon sample preparation. 

As shown in figure[9l (oc rig) is almost independent of x, including its temperature 
dependence. Considering that the volume fraction for superconducting phase (wi) is 
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Figure 8. (Color online) Relative yield of volumetric fraction for superconducting {wi) 
and magnetic (1^2) components in CaFei-xCoxAsF as a function of Co concentration 
ix). 




Figure 9. (Color online) Temperature dependence of (a) superfluid density {as = 
\/2Ss) and (b) fractional yield of magnetic phases (^2) for x = 0, 0.075 and 0.15. 
Curves in (a) are fits by models described in the text. 



nearly proportional to x, this insensitivity of Ug to x indicates that the superfluid (and 
corresponding carrier density in the normal state) stays in certain domains ("islands") 
centered at Co ions. A crude estimation suggests that the domain size may be 
ds ~ {abc/2 ■ 0.8/0.15)^/^ ~ 0.9 nm in diameter (where a, b, and c are unit cell size). In 
other words, the superfluid behaves as an incompressible fluid in CFCAF. 

Following LFAO-F, the temperature dependence of cXs in flgure[9]is compared with 
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Table 4. Physical parameters deduced from curve fitting for the data in figure [D^a) . 





Single-gap 


Double-gap 


Power law 


X 


0.075 0.15 


0.075 0.15 


X 


0.075 0.15 


a(0) ifis-') 
2Ai/A;bTc 
2A2/kBT, 

w 


20.5(8) 20.6(6) 
0.42(1) 0.44(1) 
2.4(1) 2.9(1) 


19.2(9) 20.2(5) 
0.44(1) 0.44(1) 
3.7(8) 3.9±3.1 
1.1(7) 2.3±1.5 
0.73(18) 0.5±1.3 


a(0) (/is-i) 


19.9(2) 19.8(3) 
0.45(1) 0.46(1) 
1.5(1) 1.8(1) 




0.92 1.10 


0.063 1.39 


xVNf 


0.27 1.91 



theoretical predictions for a variety of models with different order parameters. The 
weak-coupling BCS model (s-wave, single gap) again fails to reproduce the present data 
for both cases of x = 0.075 and 0.15, as they exhibit a tendency to vary with temperature 
over the region T /T^ < 0.4. A fit using two-gap model [eq. (jl])] shown by solid curves 
exhibit reasonable agreement with data with parameters listed in table HI although the 
deduced gap parameters are not necessarily consistent with the prediction of the weak- 
coupling BCS model (as the larger gap has a ratio 2Ai/A;bTc greater than the BCS 
prediction of 3.53). In such a situation, one might be concerned about the effect of flux 
pinning that often leads to a distortion of temperature dependence of a^- However, we 
found that this would not be the case for the present samples, as inferred from the fidelity 
of fs upon sweeping external field. These observations suggest that the superconducting 
order parameter is not described by a s-wave symmetry with single gap. When a power 
law, (Ts = o"o[l — (T/Tc)^], is used in fitting the data, we obtain result shown by the 
broken curves in figure [9] with respective exponents (3 ~ 1-5(1) for x = 0.075 and 1.8(1) 
for X = 0.15. As shown in the previous section, a very similar behavior is observed in 
LFAO-F near the phase boundary (x = 0.06) [27j, which is consistent with the recent 
theory of fully gapped anisotropic is-wave superconductivity [36] . 

3.3. Bai^^K^Fe2As2 

In contrast to the case of preceding two families of compounds, the substitution of 
divalent barium with monovalent potassium in Bai_2:K2^Fe2As2 leads to hole doping 
over the Fe2As2 layers. Microscopic evidence shows that the parent compound BaFe2As2 
exhibits magnetic order (SDW) below 140 K which is accompanied by a structural phase 
transition [HIES]. The situation is common to Sri_2:K^Fe2As2 [22], and it suggests that 
the electronic ground state of Fe2As2 layers in the parent compound is quite similar to 
that in LaFeAsO. However, the doping phase diagram is markedly different from LFAO- 
F, as it exhibits superconductivity over a wide range of hole content p = x/2 (per FeAs 
chemical formula) from 0.05 to 0.5 that far exceeds LFAO-F (i.e., 0.06 < x < 0.2)[13j. 
The phase diagram is also characterized by a bell-shaped variation of Tc against x, where 
the maximal ~ 38 K is attained near x ~ 0.4 {p ~ 0.2). Considering the manifold 
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nature of electronic band structure of Fe2As2 layers suggested by theoriesJUl \5U\, this 
difference may be attributed to that of the bands relevant to the doped carriers. Our 
/iSR result suggests that, while the magnetic ground state for x = 0.1 is very similar 
to that of other pristine iron pnictides, the superconducting property in the optimally- 
doped sample is characterized by gap parameters that are much greater than those of 
LFAO-F and CFCAF [29]. 

Figure [TO] shows the /xSR spectra obtained upon muon implantation to the sample 
with X = 0.1 under zero external field, where one can readily identify an oscillatory signal 
with multiple frequency components in the spectrum at 2 K. The Fourier transform of 
the spectrum indicates that there are actually two of them, one approaching to 27 MHz 
and another to 7 MHz (corresponding to internal fields of 0.2 and 0.05 T, respectively). 
This is qualitatively similar to the situation in LaFeAsO [I8l|19], and thereby suggests 
that the high frequency component corresponds to the signal from muons situated on 
the Fe2As2 layers while another coming from those located near the cation sites. The 
onset temperature for the high frequency component is close to 140 K, which is also 
consistent with earlier reports [HI [23]. 

On the other hand, it is inferred from ZF-/xSR spectra in figure [TT] that no trace 
of magnetism is found in the sample with x = 0.4. This is in marked contrast with the 
earlier report on yuSR measurements in a sample with x = 0.45 where a magnetic phase 
seems to dominate over a large volume fraction [23]. The quality of our specimen can 
be assessed by looking into bulk magnetization data, which is shown in figure [T2[(a). 
The sharp onset as well as a large Meissner fraction (> Air for ZFC) indicates excellent 
quality of the specimen. Assuming a general value for the Fermi velocity (~ 10® ki/s), 
the mean free path (Z) of carriers estimated from resistivity just above Tc is about 3.5 




0.0 0.2 0.4 0.6 0.8 1.0 2 4 6 

t/p.s 



Figure 10. (Color online) ZF-/iSR time spectra observed at 190 K and 2 K in 
Bai_2;Ka;Fe2As2 sample with x — 0.1 (Tc < 2K), where a spontaneous muon precession 
signal (mainly consisting of two frequencies) is clearly seen. Inset shows temperature 
dependence of the precession frequency. 
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Figure 11. (Color online) ZF-/iSR time spectra observed at 45 K and 2 K in 

Bai_a;Ka;Fe2As2 sample with x = 0.4 (Tc ~ 38 K). No trace of magnetic phase is 
observed. The spectra under a transverse field (TF) at 2 K is displayed on a rotating 
reference frame to extract the envelop function. Solid curves are fits using a Gaussian 
relaxation function described in the text. 



nm. Since this figure should be regarded as a lower boundary (as it is determined near 
Tc) , the actual / must be much greater than the coherence length at lower temperatures 
(^ ~ 2.1 nm, deduced from Hc2 > 75 T[51j). Thus, we can conclude that the specimen 
is in the clean limit, and that the anisotropy in the superconducting order parameter, if 
at all, should be reflected in the temperature dependence of superfluid density measured 
by yuSR. 

The /iSR time spectra under a transverse field of 0.1 T is shown in figure [TT], where 
the envelop of the damping oscillation is extracted. It exhibits depolarization towards 
zero, indicating that the entire sample falls into the flux line lattice state to exert strongly 
inhomogeneous internal field to implanted muons. The lineshape is well represented by 
a Gaussian damping, and the analysis is made by curve fit using a further simplified 
version of eq. (E]) with W2 set to zero, namely, 

1 



G^{t) = exp 



cos(27r/st + 0), (9) 



where is determined by fitting data above and subtracted from the total linewidth 
for the spectra below T^. 

The deduced linewidth, cJs, is plotted against temperature in figure [TWb). where cTg 
is the quantity proportional to the superfluid density n^. Compared with the results in 
LFAO-F [figure [3](d)] and in CFCAF [figure [Ht^a)], it is noticeable that ag rises relatively 
sharply just below Tc, and becomes mostly independent of temperature below 15 K 
(~ OATc). A curve fit using the power law, = (To[1 - {T/T^)^], yields (3 = 4.08(5), 
which is perfectly in line with the prediction of conventional BCS model for s-wave 
pairing. The gap parameter is obtained by a fit using the weak coupling BCS model 
to yield A = 8.35(6) meV and corresponding ratio 2A/ ksT^ = 5.09(4). These values 
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Figure 12. (Color online) a) Magnetization of Bai_a;Ka;Fe2As2 measured on the /iSR 
sample with x = 0.4, where data were obtained after cooling under an external field 
(FC) or zero field (ZFC). The total weight of the sample is 0.1 g, and the magnetization 
corresponding to An is obtained by using the structural parameters reported in Ref . |12j . 
b) Temperature dependence of Gaussian linwidth as = V^Ss determined by TF-/zSR 
measurement with H — 0.1 T. Curves represent the prediction of weak coupling BCS 
theory, where solid curve is a result of fit with a free gap parameter (2A/kBTc), and 
broken curve is that with 2A/kBTc fixed to 3.53 (predicted by the BCS theory). 



are consistent with the superconducting order parameter of the isotropic s-wave pairing 
with a relatively strong coupling to some bosonic excitations. 

It is inferred from recent angle-resolved photoemission spectroscopy (ARPES) study 
on BKFA with the same potassium content that the magnitude of superconducting gap 
depends on the Fermi surfaces [47\. They report Ai ~ 12 meV on small Fermi surfaces 
and A2 ~ 6 meV on the large one. We examined the consistency of our data with 
the ARPES result by employing the double-gap model shown by eq. (jlj). The curve 
fit assuming a common Tc = 38 K and a large gap fixed to 12 meV (2Ai//cbTc = 7.3) 
perfectly reproduces data in figure \T2\h) with A2 = 6.8(3) meV [2A2/A;_bTc = 4.1(2)] 
and the relative weight w = 0.30(3), where the obtained curve is virtually identical with 
that for the single gap on figure fT2](b). This, while endorsing the credibility of our data 
in terms of temperature dependence of ag, indicates that the quasiparticle excitation 
spectrum associated with multi-gap superconductivity tends to be merged into that of 
the single gap when the small gap has a large value for 2A/kBTc. 
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4. Discussion 

/iSR studies of LFAO-F have been made by a number of groups. According to tliose 
preliminary studies, no clear sign of magnetism is observed in the sample around 
X = 0.06, except for a weak relaxation observed far below Tc for x = 0.05 and 0.075 in 
ZF-/iSR spectra and an unidentified additional relaxation observed in TF-/iSR spectra 
for X = 0.075. [m [251 121] This led us to recall the sensitivity to chemical stoichiometry 
in the emergence of the spin glass-like A-phase observed near the boundary between 
the antiferromagnetic and superconducting phases in CeCu2Si2.[l5]. In addition to the 
A-phase, the present LFAO-F system exhibits a closer similarity to this classical heavy- 
fermion superconductor such as the phase diagram against pressure/doping. [36] Further 
study of the dependence of fractional yield for the magnetic phase with varying x (in 
small steps near the phase boundary) would provide more insight into the true nature 
of these phases and the mechanism of superconductivity itself that is working behind 
the coexistence/competition. 

The volumetric expansion of superconducting domains upon electron doping to 
the Fe2As2 layers by Co substitution for Fe is a remarkable feature that has no 
counterpart in high Tc cuprates. Meanwhile, this behavior to some extent reminds a 
parallelism observed in the effect of Zn substitution for Cu in the cuprates that destroys 
superconductivity: it appears that superconductivity is suppressed over a certain domain 
around the Zn atoms like a "Swiss cheese" when an extra d electron is introduced [48j. 
Despite that the effect discussed in the cuprates are completely opposite to that in 
iron pnictides, the observed "local" character of doping in CFCAF, which seems to 
come from a short coherence length C, (that probably determines the domain size, so 
that ^ ~ ds/2 ~ 0.45 nm), may provide a hint for the microscopic understanding of 
superconductivity on the Fe2As2 layers, particularly for the n type doping. 

In the meantime, the superconducting character of p type iron pnictides seems to be 
considerably different from n type ones (at least for a doping range near the boundary 
between magnetic and superconducting phases), as suggested from the behavior of 
superfiuid density observed by /iSR in Bai_j;Kj.Fe2As2. Although the temperature 
dependence of cXs in these compounds can be reproduced by the above mentioned 
double-gap model, the gap parameters shown in tables [2] and H] are considerably smaller 
than those in BKFA. The double-gap feature revealed by ARPES supports a view that 
superconductivity occurs on complex Fermi surfaces consisting of many bands (at least 
five Fe d bands) that would give rise to certain intricacy [49l [50] . Apart from the 
validity of applying the double-gap model to electron-doped iron pnictides, these figures 
suggest that the hole doping may occur in the bands different from those for electron 
doping, where the characteristic energy of the Cooper pairing may differ among those 
bands. Concerning the competition/coexistense of superconductivity and magnetism, 
there are so few data available at this stage and thus /iSR study must be extended over 
the region of x near the phase boundary (centered around p = x/2 = 0.1). Here, the 
control over the homogeneity of the specimen would be a key to proper understanding 
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of the electronic ground state. 

Finally, it would be appropriate to comment on the possible influences of impurities 
to the temperature dependence of os- First of all, they serve as the source of scattering 
for the Cooper pairs, where, depending on the corresponding length of mean free path, 
the anisotropy of the pair potential tends to be smeared out when ^ > I. This effect 
seemed to be so strong in certain high-T^ cuprates that Us might have shown a behavior 
resembling that of s-wave symmetry |52j. leading to some concern that a similar situation 
might have occurred in Bai_2^K2^Fe2As2. However, as mentioned earlier, our estimation 
on / indicates that the condition ^ -C / is satisfied in our sample, and therefore we can 
interpret the temperature dependence of cts without such ambiguity. 

Another important influence of impurities might be that they serve as pinning 
centers for flux lines. The flux pinning distorts the flux line lattice, leading to 
additional depolarization to enhance cxs. Since the pinning tends to be stronger at lower 
temperatures, it may distort CTg to exhibit stronger temperature dependence observed 
for the order parameters with nodes. (Note that it has an effect opposite to what is 
expected for the above-mentioned "scattering".) This has been tested for the case of 
CaFei_a;Coa;AsF with x = 0.075 at 2 K by shifting the external field after field cooling. 
As a result, we found that the average field in the sample indeed followed the shift 
of the external field. Thus, we confirmed that pinning was negligible at least in this 
particular case, and can conclude at this stage that the effect of impurities did not affect 
the temperature dependence of cXs seriously in CaFei_a;Coa;AsF. 

5. Summary 

It has been revealed by our /xSR experiment that superconducting and magnetic phases 
coexist in LaFeAsOi-^Fa; with x = 0.06. These two phases simultaneously develop 
just below Tc, strongly suggesting an intimate and intrinsic relationship between these 
two phases. The result of TF-/iSR measurement suggests that the superconductivity 
of LaFeAs(Oo.94Fo.o6) cannot be explained by the conventional weak-BCS model (single 
gap, s-wave). While a similar situation is suggested for the superconducting portion 
of CaFei-i^CojjAsF with x = 0.075 and 0.15, the segregation between those two phases 
are more complete, as inferred from the progressive development of superconducting 
volume fraction with increasing x. This "local" character seems to be an important 
feature common to iron pnictides (i.e., an intrinsic property of Fe2As2 layers), probably 
controlled by the short coherence length (^ 0.45 nm in CFCAF). 

Meanwhile, we have shown in a hole-doped iron pnictide, Bai_^.Ka.Fe2As2, that 
the superconducting order parameter is characterized by a strong coupling to paring 
bosons, as inferred from large gap parameters {2Ai/kBTc ^ 3.53). The temperature 
dependence of the superfiuid density, UsiT), determined by /xSR is perfectly in line with 
that predicted by the conventional BCS model with fully gapped s-wave pairing. A 
detailed analysis with phenomenological double-gap model indicates that ns{T) is also 
consistent with the presence of double-gap, although the large gap parameters make it 
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difficult to determine the multitude of energy gap solely from UsiT). 
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